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Abstract. Under proper hypotheses, we rigorously derive the Plateau angle 
conditions at triple junctions of diffused interfaces in three dimensions, starting 
from the vector-valued AUen-Cahn equation with a triple-well potential. Our 
derivation is based on an application of the divergence theorem using the 
divergence- free form of the equation via an associated stress tensor. 



1. Introduction 

We consider the problem of determining contact angle conditions at triple junc- 
tions of diffused interfaces in three-dimensional space via the elliptic vector-valued 
AUen-Cahn equation 

(1) Au - V„VK(u) = 0, 

for maps it : — > and a triple-well potential : R^ — > M. Equation ^ is the 
vector analog of the well-known scalar elliptic equation 

(2) Au-W'{u) = 0, 

for It : M'' R and a bistable potential : R — > R with two minima, which was 
introduced by Allen and Cahn ^ in the context of antiphase boundary motion. 
Here, u is an order parameter that denotes the coexisting phases of the phenom- 
enon of phase separation. The vector-valued version was considered by Bronsard 
and Reitich |14j as a generalization for more than two phases (see also Rubinstein, 
Sternberg, and Keller [25]). Note that both equations are elliptic versions of corre- 
sponding evolution problems that involve a small parameter e, which denotes the 
thickness of interfaces, and that they are Euler-Langrange equations of energy func- 
tionals, whose minimizers are related to minimal surfaces (see Modica and Mortola 
[111, Modica [21], and Baldo [11] )• 

For the problem of contact angles in the case of soap films in three dimensions 
the classical Plateau angle conditions state that 

(1) three soap films meet smoothly at angles of 120 degrees along a curve, 

(2) four such curves meet smoothly at angles of about 109 degrees at a point. 
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The above laws hold in the isotropic case of soap films, which corresponds to systems 
of minimal surfaces (cf. Dierkes et al. [HI §4.15.7]). The angle of 120 degrees in 
^ is the angle whose cosine is — ^, which is exactly the angle in isotropic triple 
junctions in two dimensions, while the angle of about 109 degrees in ([2]) is the 
angle whose cosine is — ^ (the so-called Maraldi angle). In the anisotropic case of 
mixtures of immiscible fluids, the angles above are not always equal and depend 
on the surface tension coefficients of each fiuid, as in systems of constant mean 
curvature surfaces. In this case, the angles at the four triods determine the angles 
at the point singularity. 

Going back to equation ([ij, in [M] the authors proved short-time existence and 
uniqueness of solutions to the parabolic vector- valued £-problem in the case of triple 
junctions on the plane and, using formal asymptotics, showed that the solution 
partitions its domain in regions where it is approximately equal to one of the three 
minima of the potential. Then, by blowing up around the triple junction and 
at a slow timescale they also showed that, to leading order, the solution is time- 
independent and that it solves the stationary equation ([1]). Moreover, near the 
junction the interfaces are fiattencd out and from the matching conditions of the 
asymptotic analysis it follows that the limit along directions parallel to an interface 
depends only on the distance to the interface. 

From the rigorous viewpoint, in two dimensions, we also mention the existence 
results of Bronsard, Gui, and Schatzman jl3j for triple- well potentials and of Alama, 
Bronsard, and Gui [T] for potentials with multiple minima, both under assumptions 
of symmetry (see also Saez Trumper |26|). In three dimensions, Gui and Schatz- 
man |19| proved existence for potentials with four minima, again under symmetry 
assumptions, while Alikakos and Fusco [7J [31 [T7] proved existence in arbitrary di- 
mensions for multiple-well potentials that are invariant with respect to a finite re- 
flection group of symmetries. Such solutions partition the domain space in regions 
where they arc approximately equal to the minima of the potential, separated by 
flat diffused interfaces. Concerning the partitioning problem in three dimensions, 
Taylor |28| classified all possible partitioning cone configurations and showed that 
the only locally minimizing ones are a single flat surface, a triod of flat surfaces 
meeting along a line at 120 degree angles, and a complex of six flat surfaces meet- 
ing with tetrahedral symmetry at a point, as in Plateau's laws. To our knowledge, 
complete rigorous results are not available in arbitrary dimensions and we refer to 
Morgan [231 Ch. 13] for a review. 

For the rest of this note we restrict to the three-dimensional case and consider 
uniformly bounded entire solutions to ([T]) , as constructed in [71 [Sj [17] . Taking the 
potential W to have three minima, such solutions partition M.^ in three regions that 
are separated by three interfaces which intersect along a line that we call the spine 
of the triod. In the symmetric case, this corresponds to one of the two singular 
minimizing cones for the associated Plateau problem. 

For the problem of determining contact angles, another result in [IJ (and also 
in Gui [18]) is the derivation of the law 

sin 01 _ sin 02 _ sin 03 

tT23 Cai (Tl2 

for the angles 0i, 02, 03 between interfaces at a planar triple junction, where the 
a's are surface tension coefficients at the interfaces between neighboring phases. 
The authors have called this relation Young's law, after Young [53], who studied 
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the contact angle of a droplet on a solid flat substrate. Another name is Neumann's 
triangle, after Neumann |24) . who used this relation for triple junctions of fluidic 
interfaces (cf. Minkowski [501 P- 567]). In three dimensions, the angles at a triod 
also obey Young's law, while the angles at a quadruple junction are a geometric 
consequence of the angle conditions at the four triods that form it (see Bronsard, 
Garcke, and Stoth |12j for the calculation), that is, given a quadruple-junction 
configuration. Plateau's second law follows from the first, which is the isotropic 
version of Young's law. 

Our goal in the following is to rigorously derive Young's law for triods of interfaces 
in three dimensions as a property of solutions satisfying the next two hypotheses. 
(These will be made precise in Section[2]) We note that in the context of symmetric 
solutions, the hypotheses below are theorems (see [71 [3l \YJ\ [8]). 

Hypothesis 1. In the interior of regions and along directions extending from the 
spine to infinity, solutions converge exponentially to the corresponding minima of 
the potential. 

Hypothesis 2. In the interior of regions and along directions extending to infinity 
while being parallel to an interface, solutions converge pointwise to connections, 
that is, maps with argument the distance to the interface and with the property of 
connecting the minima of the potential at plus and minus infinity. 

The derivation makes use of the fact that equation (IlJ is a divergence-free con- 
dition for a certain stress tensor, which appeared in [5] in this context (see also [S] 
for further comments on its origins). For the planar analog, the derivation in jl4j is 
based on formal asymptotics, while the derivation in jl8[ is closer to our spirit but 
uses Pohozaev-like identities instead of the stress tensor (see [5] for the connection 
between the two). For the three-dimensional problem, related results appear in the 
theses [TU] and [H]. 

The rest of the present note consists of two sections. In Section |21 we set up 
the problem, stating all necessary assumptions, and formulate it in divergence-free 
form via the stress tensor. In Section |3l we apply the divergence theorem for the 
stress tensor on a ball, which is then blown up after a proper slicing. This slicing 
involves a surgery around the singularity which appears at the intersection of the 
surface of integration with the spine and a breaking up of the remaining part in 
a way that utilizes the hypotheses on the solutions at infinity, that is. Hypothesis 
in far from interfaces and Hypothesis [21 at fixed distances from them. This yields 
Young's law in the form of a balance of forces relation for the conormals of the 
three interfaces, which is equivalent to ([3]). 

2. Statement of the problem and preliminaries 
We start with equation ([1]) 

Au - VuW{u) ^ 0, 

for u : K3 ^ R3 g^nd : R3 ^ R, where WuW{u) = {dW/dui,dW/du2,dW/du3). 
The potential W is taken to be of class C^, nonnegative. and with three nondegen- 
erate global minima at points ai, 02, 03, that is, W{ai) — 1^(02) — W{a3) — 0, 
with W{u) > otherwise. Moreover, we ask that W satisfies a certain coercivity 
assumption, that is, that there exists M > such that W{su) > W{u), for s > 1 
and \u\ = M. Finally, note that we do not make any symmetry assumptions on W. 



4 



N. D. ALIKAKOS, P. ANTONOPOULOS, AND A. DAMIALIS 



As explained in Section [U we consider solutions that partition the domain space 
in three regions via a triod of diffused flat interfaces. We distinguish three regions 
Ci in M'^, for i = 1, 2, 3, such that the region Ci contains the minimum a^, with Tij 
being the interface that separates Ci and Cj (with Tij = Tji). For each region Ci 
we have that if x e C;, then Ax e C;, for A > (cone property). 

We choose coordinates as follows. We take the origin on the spine, which we 
identify with the X3 axis, and we further identify interface T12 with the half plane 
xi = 0, X2 > 0, X3 £ E, such that xi is the distance to ri2. We also recall here the 
spherical coordinates in three dimensions, that is, 

Xi = r cos 01 sin 02, X2 = r sin^i sin02, 2:3 = r cos 02, 

for an azimuthal angle 61 £ [0, 27r), a polar angle 6*2 G [0,7r], and for r > 0. In 
terms of the azimuthal angle 61, the interface ri2 lies at = ^. 
The uniformly bounded entire solutions we consider satisfy 

(4) \uix)\ < C, 

globally in R'^. Using this bound and linear elliptic theory, we also have the uniform 
bound 

(5) |Vu(x)| < C, 

again globally in M^. 

For such solutions we have two hypotheses. The first one concerns the fact that 
solutions converge exponentially to the corresponding equilibrium in the interior of 
each region. This has been verified under assumptions of symmetry on the potential 
by several authors (see [131 [191 [71 [Sj [17] ) and we postulate that it holds for general 
potentials. 

Hypothesis 1 (Exponential estimate) . In the interior of the region Ci there holds 
that 

(6) |w-a,| <e-'*''^*(^'^^'', 
where dCi = Ui^jTij 

(We use the notation X <Y ioi the estimate X < CY, where C is an absolute 
constant.) 

The second hypothesis is that along directions parallel to interfaces solutions con- 
verge to one-dimensional heteroclinic connections Uij , which connect the equilibria 
tti, aj at infinity, in the sense that 

lim Uij{r/) ^ Qi and lim UiJri)=aj, 

where ri is the distance to the interface Tij. We refer to [27l [6l [4] for further 
information. Hypothesis [2] has also been verified for symmetric potentials under 
the additional assumption of uniqueness and hyperbolicity of connections [8] . 

Hypothesis 2 (Connection hypothesis). Along directions parallel to an interface 
Tij solutions converge pointwise to a one- dimensional connection Uij^rj) with argu- 
ment the distance to the interface, that is, 

(7) lim u{x) ^ Uij , for fixed 1] := dist{x, Tij). 
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These limiting functions are solutions to the associated Hamiltonian ODE system 

with the property of connecting the minima of W at infinity. We define the action 
of such a connection to be the nonncgative quantity 

(8) = ^(f/,,) Qlt/^jf + WiU,,)^ d77, 
and also note that connections satisfy the equipartition relation 

(9) IM = W{U.,,). 

We will now reformulate equation (HI) via its associated stress tensor (see [5] for 
more information). We define the stress tensor T as 

(10) T,j{u) = ■ u^, - % + ^("^ 

for maps u : M" — > M™, where u^i ~ du/dxi and where the dot denotes the Eu- 
clidean inner product in M™. In three dimensions (that is, for n = 3) it is a 3 x 3 
symmetric matrix 

|2 |„, 12 L, 12 




W k,ir-|".2r-|M.3| -2VF(w) 2^.1- M.2 2w.i-",3 

T(u)=- 2W.2-U.1 |M2p-|ui|'-|".3|'-2iy(u) 2W2-W3 



with the property 

(11) divT = (Vu)^(Ait- V„M^(m)), 

that is, T is divergence- free when applied to solutions of equation ([1]). 

We also note that T is invariant under rotations of the coordinate system, that 
is, it transforms as a tensorial quantity. To see this, consider an orthogonal trans- 
formation Q and a new coordinate system x' = Qx. Letting u' be the map acting 
on the new coordinates, with u'[x') = u(x), the chain rule gives that its gradient 
is transformed via V'u' = QVu, where the prime denotes that the derivatives are 
taken with respect to the new coordinate system. Then, for the transformed tensor 
T', which is given by the similarity transformation 

r = QTQ^, 

due to the form of the components in (fTO|) and the continuity of W there holds 

Triu')^u[,-u[^-d.,(^l\Vuf + Wiu': 

where again the prime denotes that the tensor is calculated in the new coordinate 
system. That is, the transformed tensor has exactly the same expression as the 
original one, except for the fact that it acts in the new coordinate system. 

Finally, we state without proof two lemmas that will be used in the following. 
The first is a consequence of Hypothesis[T]and linear elliptic theory, while the second 
follows from Hypothesis [5] and the Arzela-Ascoli theorem. 

Lemma 1. Solutions of equation ([T]) satisfy the gradient estimate 

(12) |Vu(a;)| <e-'^''''("'^^'\ for x e C 
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Moreover, a similar estimate holds for the potential W{u), that is. 



(13) \Wiu{x))\ < e-di^'t^'^c^-), for x € Q. 

Lemma 2. For solutions of equation ([T]) the following pointwise limits hold. 

(14) limit. i(a;) — U{xi), as X2 — > +oo, x^ — > +oo, 

(15) lunuj{x) = 0, as X2 — > +oo, x^ — > +oo, for z = 2, 3, 



where without loss of generalit-^ we considered a coordinate system such that xi is 
the distance to an interface. 

3. Derivation of Young's law 

In this section we will prove the following theorem. 

Theorem. For the contact angles at the spine of a triod of intersecting interfaces 
Ti2, T23, Fai, Young's law holds in the form of a balance of forces relation, that is, 

(16) (Tl2iyi2 + CF23V23 + CT^lV^l = 0, 

where atj is the action of the connection Uij of each interface Tij and Vij the 
corresponding unit conormal, that is, a unit vector that is tangent to Tij and normal 
to the spine. 

Proof. Since the solutions of equation ([T]) which we consider are constructed as 
minimizcrs over balls (see [71 [3j I17| for the variational setup of the problem) , we 
take a ball Bn centered at (0, 0, 2R) in order to apply the divergence theorem on it 
using ([Tl]) . that is, the fact that the stress tensor T is divergence-free. This gives 

(17) 0=4/ diyTdx^^ [ TudS, 

R J Br R JoBr 

where v is the outer unit normal to the boundary dBu. In what follows we will 
study the limit 

1 f 

lim — / TvdS 

in order to utilize the hypotheses on the solutions at infinity. Note that we chose 
the center of Bj^ in such a way so that for {xi, X2, x^) G dBj^, we have 2:3 7^ and 
X3 — >■ +00 as R ^ +(X). 

The complication in applying the divergence theorem in our problem is that 
the surface of integration intersects with the spine at two points, where we have 
no information on the behavior of solutions. In our setup these arc the two poles 
of Bn, at (0,0, i?) and (0,0, 3i?). To circumvent this, we perform a surgery by 
choosing two appropriately sized spherical caps around the poles. To this end, let 
-02 (i?) be a small polar angle that defines the spherical caps (see Figure[T]) for which 
we require that there holds 

(18) i? sin ?/'2 (i?) +00, as i? ^ +00, 

such that the distance of the boundary of the cap to the spine grows as i? — >■ +00, 
which also yields that the geodesic radius Rip2{R) of the cap grows as i? — +00. 
Moreover, we require that 

(19) Rip2{Rf ^0, as i? ^ +00, 

^Duo to the invariancc of the Laplacian under rotations and the continuity of W and V^W. 
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Figure 1. The sphere Bn centered at the spine, with two caps of 
polar angle tp2 {R) and a strip at geodesic distance Rtjji (R) around 
the intersection of the interface with Bj^. 



such that the renormalized area of the cap (in the sense below) shrinks as i? — > +00. 
To see this, first note that condition also yields that 



(20) 



ip2{R) -> 0, as i? -> +00. 



The renormalized area of a cap can be easily calculated as a surface integral using 
spherical coordinates, that is, 

sin 6*2 de'2 de*! = 27ri?( 1 - cos V'2 ) ■ 



1 
R 



dS 



1 

R 



Using (gni) we have that (1 - cosi'2{R)) = 0(?/'2(i?)^), which via ^ gives 
(21) ■• ^ 



lim - / dS* = 0. 



To sum up, we choose the size of a cap to be small enough so as not to matter in 
the integration, but at the same time large enough so we always stay away from 
the singularity. Then, for the integral of Tv on such a cap we have the estimate 











j TudS 


<!/ 


R 


' cap 


^ J cap 



1 

'r 



d5. 



cap 



where we used the bounds (HJ, ((5)), and estimate of Lemma [T] for bounding |r| 
by a constant. Using now we finally have 



lim — 
R 



TvdS = 0. 



For the remaining part of the sphere, we will work separately for each interface 



that intersects it. For the interface ri2, which lies at azimuthal angle 9i 
work with the slice 



we 



^={(^1,^2,0 



< 
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for a fixed angle 5, such that no other interface intersects with the slice. To study 
the limit 

lim — / Tv dS, 



fi^+oo R Js\(5ncaps) 

wc distinguish two parts in S \ {S D caps), a neighborhood around the meridian 
at the intersection of the interface with the sphere and the rest. We take the set 
J\f C S \ {S O caps), such that J\f is the strip that is contained between two planes 
parallel to ri2, one in Ci and one in C2, and at equal distance Rsin'ijji{R) from it 
(see Figure [1]). For the angle ipi{R) we require that it satisfies 

(22) t/-! (i?) < {R) , with V2 sin V-i (i?) < sin i/'2 {R) , 

such that the azimuthal angle -01 that defines the width of the strip JV at the 
equator of Bji is strictly smaller than the polar angle 'ip2 that defines the caps. 
This condition forces A/" to be a subset of S' \ (S* n caps). Moreover, we require that 

(23) i?sin7/'i(i?) -> +00, asi?— >+o3, 

such that the distance of the interface to the two planes that define J\f grows as 
R — > +00. From conditions (|22p and (PU)) we also have that 

(24) iPi{R) -> 0, as i?^ +00. 

An example of angles ipi , V'2 that satisfy all the above requirements is 

MR) ^ R^'^''" and i}2{R) = 

with condition ([^^ holding true for R > 1025 for this particular choice. 
Given the following decomposition of the set S\{S O caps), 

S\{Sn caps) =Af U ((5 \ (5 n caps)) \ N) , 

we have the estimate 
1 
R 



TvdS 

(S\(Sncaps))\A^ 



< 4 / \T\\v\dS 



R 



(5\(Sncaps))\A/' 



(25) <- [ g-flsi„0i(_R) 



R 



(S\(Sncaps))\A^ 



using estimates p^ . ([T^ from Lemma [T] for estimating |r| by the exponential 
and since the domain of integration is of order 0{R^). Finally, taking the limit as 
R +00 and using condition (p3)) . we have that 



lim 4 / Ti^dS^O. 

R J(5\(5ncaps))\AA 

We turn now to the last part, which is the integral on the strip M. We parame- 
trize Af as the graph of 



Mxi,X3) = ^i?2-x2-(x3-2i?)2. 



for xi e (— i?sin?/)i, _Rsin V'l), X2 G {R\/ sin^ — sin^ ipi,R), X3 e {2R~ Rcos^2, 
2R + RCOS1P2), where xi = 2:2 = = is the origin, which is not the center of the 
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sphere. We set yi = xi, y2 = X2, and 2/3 = — 2R, such that (2/1,2/2,2/3) G dBa 
with yl + yl + y3 ~ R^- Then, Af is the graph of 



/ii(yi,2/3) = \Jr^ -yl-yl 

for 

2/1 G (— i?sin-0i, i?sin-0i) 



2/2 e sin^ i'2 - siv? ipi, R) 
ys e (-i?cos'(/;2, -Rcos-02) 
For the surface element we calculate 



(26) 



dS=^l+y-^j + (^"^ j d2/3d2/i = — d2/3d2/i, 
where 2/2 = ^/i?'^ — 2/1 ^ 2/3 > while the outer unit normal is 

y 

R 



dfnV , fdfR\\ , R 



Using this parametrization, the integral on J\f is written as 

1 /• 1 (--R sin -01 i-R COS 412 

(27) - TvdS^- / T(«)^d2/3d2/i, 

-tt JA/" ^ J-_Rsin i/Ji J-_Rcosi/J2 ^2 

where 

u(xi,a;2,a;3) = w(2/i, 2/2, 2/3 + 2i?) =; w (2/1, 2/2, 2/3), 
and i',i(2/) = w_i(a;), for i = 1, 2, 3, so T{v{y)) = T{u{x)). 

To take the limit as i? — > +00 in equation (|27p and apply Hypotheses [U [21 
we use Lebesgue's dominated convergence theorem. The components of the vector 
quantity to be integrated on the right-hand side of ([27]) are given by 



(28) (t{v)^] =T,,(i;)^, for ^ = 1,2,3, 

V 2/2/, 2/2 

using the summation convention. To check whether dominated convergence applies 
for each component, we write the corresponding integral as 

^X[-flsinVn,i?sin0i] / T^j{v)—dy3 dyi, 

3 V-n- J-Rcosi,2 2/2 / 

where X is the characteristic function, and we would like to show that the quantity 
in parentheses is dominated by some integrable function. Using the estimates 
and ^ of Lemma[Il in TV there holds \Tij{v)\ < e^l^'il, which gives the estimate 



1 f'''°'^' rr. , ,yj 



nX[-i?sin^/>i,_Rsini/Ji] / Tij{v) dy3 

-n- J-Rcosi'2 2/2 



1 /--R cos 1^2 I I 

<^/ m,(z;)i^d2/3 

^ J-Rcos4>2 2/2 

<1 r^™^'^^ 



Setting 



RJ-Rcosip2 2/2 



e-i«ii^d2/3 

^ J-Rcosi>2 2/2 

\^J-flcos02 2/2 



d2/3, for i = 1,2,3, 
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we argue 



wc aiguc that Ij < 2, for j = 1, 2, 3, and, as a consequence, dominated convergence 
applies in the hmit R — > +oo. 

For J = 1, using the extremum values of j/i, y2 in the intervals (|26p and condition 
(l22l) . we have 



1 



d2/3 < 



/•i?, cos 'i/;2 



i? sin -01 



1 /•«™'^'/'2 gjj^^^ 1 I-B.C0S4,2 

< — I -. — - = p / dya = 

^ J —R cos 1p2 



^ J -RcosiIj2 R\/ sin^ ?/'2 — sin^ i/'i 
2 cos ^^2 < 2 



"02 — sin^ 01 from ((22 



since sin < \/sii?^ 
For .7 = 2, we have 

h^^ —dyz^- d2/3 = 2cosV2<2, 

^ J-Rcosip2 2/2 ^J-flcos i/'2 

since y2 > in A/" 



ice 2/2 > in A/". 
Finally, for .7=3 we change variables to 



2/3 = ^2^3, with dy3 = ^ dys, 
K It 

and, using the extremum values of yi from (|26p . we estimate 



1 



^3 = 15 



,i?COS^2 ly^l 

dys = 



< 



.COSV.2 1^3 



^cos^2 ^1^31 



-cosi/;2 2/2 J -cos -4:2 y/R'^ - yl - R'^y'l 



:dy3 



i/3 = 2 



-cos,/;2 ^1 _ sin^ _ y2 



dy3 



/•cos -02 



ys 



1 — sin^ 01 — y| 



dy3, 



since the function 
integral to get 

^COS ^/l2 

2 / 



I y3 1 / y^l — sin^ "01^2/1 is even. We explicitly calculate the last 

/■COS tl)2 



1 - sin^ V>i - y| 



dy3 = -2 



1 - sin^ i'l - y| ) dy3 



= — 2y^ 1 — sin^ tAi — cos^ 7/>2 + 
= 2 cos 0^1 



2 J 1 - sin^ -01 



— 2\J sin^ 012 — sin^ 0i < 2 cos0i < 2. 

in ([27]) . wc distinguish 



To conclude with the calculation of the limit as i? ^ +00 in ([27]) . wc dif 
the following limits in J\f, as consequences of Hypothesis [5] and Lemma [H 

limfl_>+oo v{y) = limx2^+oo u{x) = Ui2{xi), 

limn^+oo WaCy) = limx2->+oo u.ii^) = t^i2(a;i), 

limH->+oo w,2(y) = limx2->+oo ",2(2;) = 0, 

limH->+oo 'y,3(y) = lima;2->+oo ",3(2;) = 0. 

Using the extremum values of the intervals in ((26| and condition (p2)) , we also have 



(29) 



that in A/" 



i? sin -01 



^ i^i/sin^ 0^2 - sin^ 0i ^ i?sin0 



1 sin -01 -r^ 
— -. — ■ > 0, as i? +00, 
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1 RCOS'4'2 1 C0S?/'2 „ „ 

< , , = < -> 0, as i? ^ +00, 

R i? Vsii7^/^2^^siiiVL ^ ^'1 



where in the last limit we also used condition (j23p for the limit of the denominator. 
Using now ([50)1 . ([5T|) and since the elements of T are bounded by a constant from 
dH) , (O , and estimate of Lemma [TJ we have 

lim ^r(w)— = lim ^{Ti2,T22,T32)^ , 

R^ + oo H 1/2 R^+OD H 

that is, only the components for j = 2 in (j28p do not vanish in the limit. But, using 
(l29l). we further have that 



lim — 2^12= lim — vi-V2=0 and lim — 732= lim — V3-V2=0- 

R^+oo R R^+OD R ' ' R^+oo R ' R^ + oa R ' 

Finally, 

fj?sim/)i pR cos 11^2 



\im - Ti.dS=i Jim - / T22(«)d2/3dyi (0,1,0) 

Plugging in the component r22 into the last integral, we calculate the limit 

-1 /'i?sin'0i />i?. cos?/^2 1 

1 / / ^ ' ' ' |2 I |2 



/•i? sin -01 


pRcos lp2 




J sin ^ 


' —R, cos '(/J2 





T 



^hm^ - / / _ ( 1^ 2 1' - |« 1 1' - It'.a I' - 2iy (i;) ) d2/3 dyi 

via the change of variables = Ry^ , which gives 

/Rsinipi (.cosi/;2 l , i . 

/ 7^[VK2? - \vs? - -^k,3p - 2W{v)) dys dyi, 

-flsin V'i J- cos 02 ^ ^ ^ ^ 

with a slight abuse of notation for w.3. Passing the limit inside the last integral and 
using conditions (PH]) and ([25]) . the limits in give 

1 a 



^- (^2|C^i2(yi)r + Vl^(C/i2(2/i))j d^3cl2;i 
= -2 j Q|[/i2(yi)p + W^(C/i2(yi))) dyi 

= -2(712. 

Thus, we have shown that for the (5-slice S around the interface ri2 there holds 

lim 4 / = -2(Ti2i'i2, 



R J g 

where 1^12 = (0,1,0)^. 

Since the stress tensor T is invariant under rotations, we can apply the same 
procedure for the other two interfaces for appropriately rotated coordinate systems 
and appropriate 5-slices (in order to cover the whole sphere) to get 

0'12I^12 + 0-23i^23 + 0-311^31 = 0, 

using (jl7p . where the Vij^s are the conormals of the corresponding interfaces F^. 
This concludes the proof. □ 

We remark that the balance of forces relation (|16p is equivalent to Young's law 
(131). This can be easily deduced by multiplying (|16p with the unit normal of each 
interface. 
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